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Abstract 

We compute the Form Factors of the relevant scaling operators in a class of inte- 
grable models without internal symmetries by exploiting their cluster properties. 
Their identification is established by computing the corresponding anomalous di- 
mensions by means of Delfino-Simonetti-Cardy sum-rule and further confirmed 
by comparing some universal ratios of the nearby non-integrable quantum field 
theories with their independent numerical determination. 



1 Introduction 



In this work we present a detailed investigation of the matrix elements 



(A,...,A0 = (0\O(0)\A ai (p 1 )...A an ((3 n )) 



(1.1) 



(the so-called Form Factors (FF)) in a class of integrable two-dimensional quantum 
field theories. Our specific aim is to check some new theoretical ideas which concern 
the relationships between three different regimes which two-dimensional quantum field 
theories may have, namely the ones ruled by conformal invariance, integrable or non- 
integrable dynamics. 

Conformal Field Theories (CFT) and the associated off-critical Integrable Models 
(IM) have been extensively studied in the last years: as a result of these analyses a good 
deal of information has been obtained particularly on correlation functions of a large 
number of statistical mechanical models in their scaling limit and on physical quantities 
related to them (see for instance [1-8]). In this context, a crucial problem often consists in 
the determination of the spectrum of the scaling operators away from criticality, namely 
their correct identification by means of the set of their Form Factors. This is one of the 
issues addressed in this work. 

Form Factors also play a crucial role in estimating non-integrable effects. Let us 
first recall that the above CFT and IM regimes cannot obviously exhaust all possible 
behaviours that statistical models and quantum field theories can have since typically they 
do not possess an infinite number of conservation laws. This means that in general we have 
to face all kinds of phenomena and complications associated to Non-Integrable Models 
(NIM). The scattering amplitudes and the matrix elements of the quantum operators will 
have in these cases a pattern of analytic singularities due both to the presence of higher 
thresholds and to the appearance of resonances. A first step forward in their analysis has 
been recently taken in ref. |J where it has been shown that some interesting examples of 
NIM may be obtained as deformations of integrable models. The action of such theories 
can correspondingly be written as 




(1.2) 



1 



Aint being the action of the integrable model. Since the exact expressions ( |1.1|) of the 
Form Factors of the integrable theories are all assumed calculable, in particular the ones 
of the fields *&i(x) entering eq. (|1.2|), one is inclined to study the non-integrable effects by 
using the Born series based on the Form Factors. Although at first sight this still remain a 
difficult task (and generally, it is indeed so), there may be favorable circumstances where 
the analysis simplifies considerably. For instance, as far as there is only a soft breaking 
of integrability, it has been shown in || that the complications of the higher terms in the 
series can often be avoided since the most relevant corrections only come from the lowest 
approximation. If this is the case, one can extract an important amount of information 
with relatively little effort: a significant set of physical quantities to look at is provided 
for instance by universal ratios, like the ones relative to the variations of the masses or 
of the vacuum energy density S vac : if the breaking of integrability is realized by means 
of a single field \l/(x), those are expressed by 

5m t _ mf F£(m) 



5m, „£ 0) 

(1.3) 

SS vac (0 | g | 0) 

where mf refers to the (unperturbed) mass spectrum of the original integrable theory. 
It is thus evident that also to estimate the non-integrable effects associated to a given 
operator ty(x) one must face the problem of correctly identifying its FF's. 

Two new results on the relationship between CFT and IM have been recently derived 
by Delfino, Simonetti and Cardy JTU]. Briefly stated, the first result consists in a new 
sum-rule which relates the conformal dimension A^ of the operator <p(x) to the off-critical 
(connected) correlator (6(x)0(O)) c , where 9(x) is the trace of the stress-energy tensorQ 

A* = -— L_ [ c / 2 x(9(x)0(O)) c . (1.4) 



4tt(0) 

This sum-rule is closely related to the analogous expression for the conformal central 



1 The sum-rule in the form of eq. (|l.4[ ) may be violated by effect of renormalization of the operators 
outside the critical point, as clarified in the original reference JLQ]. This is however not the case for the 
field theories and the operators considered in this work. 
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charge c [11 



— / d 2 x | x | 2 (e(ar)6(0)) c . (1.5) 
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Equations (|1.4|) and ( |1.5| ) express elegant relationships between conformal and off-critical 
data, but more importantly, they provide very concrete and efficient tools to characterise 
the scaling limit of the off-critical models. 

As for the second result, it has been suggested by the aforementioned authors of ref. 
IIJ , that the Form Factors of the relevant scaling fieldsQ of an integrable field theory - in 



absence of internal symmetries - are in one-to-one correspondence with the independent 
solutions of the so called cluster equations 



hm Fl a2 _ ak<ak+i _ ak+i pi, e 2 , . . . , e k , a + e k+1 , . . . , a + e k+l ) = (1.6) 

i 



/£V - al,a 2 ,...,a k (#1, #2, • • • , 9 k ) ^ afc+1 . .,a fe+; (0fc+l> • • • , Ok+l) 

These equations can be imposed on the Form Factors in addition to the customary func- 
tional and residue equations which they satisfy (see in this respect also |], If this 
cluster hypothesis is valid, we would have a clear method to identify the matrix elements 
of all the relevant operators, at least in the case of theories without symmetries. It must 
be stressed that until now this task has often been a matter of keen guess-work and 
mostly based on physical intuition. 

It turns out that a check of the above cluster hypothesis provides a well-suited forum 
for testing several theoretical aspects. In fact, the most direct way of confirming the 
above idea is firstly to solve the general functional equations of the Form Factors with 
the additional constraints of the cluster equations (|1.6|) and to see whether the number 
of independent solutions equals the number of relevant fields in the corresponding Kac 
table. If the above check turns out to be positive, one may use the sum-rule ( |1.4| ) in 
order to achieve the correct identification of the (supposed) primary relevant operators 
<f>i'. from the values of the partial sums one can in fact infer the value of the anomalous 
dimension and correspondingly recognize the operator. An additional confirmation may 
also come from the employment of eqs. (|1.3f) relative to non-integrable field theories. In 



2 Hereafter we are using the short expression "scaling fields" to actually denote the off-critical operators 
which reduce to the scaling fields in the conformal limit. 
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fact, one can regard the primary field <pi{x) under investigation as that operator which 
spoils the integrability of the original theory and therefore compare the predictions ( |1.3|) 
based on its Form Factors with their independent numerical determinations which may 
be obtained by means of the truncation method ||13|1 . Note that a successful test of this 
kind could also be interpreted the other way around, namely as a further proof of the 
effectiveness of the formulas ( |1.3| ) in estimating non-integrable effects. 

The models on which we have chosen to test the above considerations are integrable 
deformations of the first representatives of the non-unitary conformal series^ -M(2, 2n + 
1), n > 2. They belong to the class of universality of solvable RSOS lattice models d 



la Andrews-Baxter-Forrester although with negative Boltzmann weigths [14], |lq| : their 
simplest example is given by the quantum field theory associated to the so-called Yang- 
Lee model which describes the distribution of zeros in the grand canonical partition 



function of the Ising model in a complex magnetic field [[U| [17|]. These models do not 
have any internal symmetry and all their fields are relevant operators: hence, they are 
ideal for our purposes. Moreover, the nature of their massive and conformal phases is 
simple enough. The price to pay for their relative simplicity is however the presence of 
typical non-unitary phenomena, as imaginary coupling constants or negative values of 
the anomalous dimensions and central charge, together with the anomalous poles in the 
S-matrix which induce an unusual analytic structure in the Form Factors 0, |17|, |I8| . 

The paper is organized as follows. In Section 2 we discuss the general strategy which 
can be employed in order to compute the FF's of the relevant operators in the integrable 
deformations of the models -A4(2, 2 n+1). In Section 3 and 4 we present a detailed analysis 
of the FF's of the models -A4(2, 7) and -M(2, 9), which are the first non-trivial examples 
on which to check all the theoretical ideas discussed above. In fact, for the first model 
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The conformal weigths and central charge are given respectively by 

(o-l)(2n-a) 

Ai a = ; ; , a = 1, 2, . . . , in 

2(2n + l) 

2(6n'-7n + l) 
2n + 1 
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A4(2,5) the cluster hypothesis is easily verified: the only solution of the Form Factor 
equations is the sequence of functions determined in ref. |7j which indeed fulfill the cluster 
equations ( |l.(j| ) and are easily identified with the matrix elements of the only relevant field 
of the Yang-Lee model. The two models Ai(2,7) and Ai(2,9) represent, somehow, the 
best playground for our purposes because they give rise to integrable models under each 
of their possible (individual) deformations and also because they optimize the size of the 
lenghty numerical output which we present for the solutions of the non-linear equations. 
Moreover, although there is in principle no obstacle to extend the analysis to all the 
models Ai(2,2n + 1), these are the simplest cases from a computational point of view 
since the larger is the value of the index n the higher is the order of the system of algebraic 
equations to be solved for determining the Form Factors. Finally, our conclusions are 
in Section 5. Two appendices complete the paper: Appendix A gathers all important 
formulas relative to the parameterization of the two-particle Form Factors and Appendix 
B collects the S-matrices of the models analysed. 



2 Outline of Our Strategy 

In this section we discuss the general strategy needed in order to obtain the Form Factors 
of the scaling primary fields of the integrable deformations <pi ^ of the conformal models 
A4(2,2n + 1) (hereafter denoted by the shorthand notation [M(2,2n + l)](i,&)). The 
deforming field <p\^ can be one of the operators 01,2, 0i,3 or possibly some other primary 
field which gives rise to an integrable deformation. 

The starting point in the computation of the Form Factors is the correct parameteriza- 
tion of the two-particle ones which is given detailed in Appendix A. This is a non-trivial 
task in the case of non-unitary models for the reason that the exact S'-matrices of these 
models are usually plagued by a pletora of anomalous poles |T8[. By this we mean for 



example simple poles which are not related to any bound state, or, more generally, any 
poles which do not have apparently the standard diagrammatic interpretation of refs. [[IS 



Consider for example the S'-matrices listed in the tables of Appendix B relative to the 
integrable deformations of the models M.(2, 7) and M(2,9) where the anomalous poles 
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have been labelled with £>, T> or *. The origin of these poles may be explained according 



to the ideas put forward in |20|]. In particular, poles of type B and T> are due to multi- 
particle processes of the kind described respectively by the "butterfly" and "dragonfly" 
diagrams drawn in Figures 2 and 3 respectively. These multi-loop processes induce in the 
S'-matrix simple poles rather than higher order ones because the internal lines of these 
diagrams cross at relative rapidity values relative to some zeros of their corresponding 
two-particle S-matrix element: this gives rise to a partial cancellation of the poles. 

The adopted parameterization for two-particle FF's is directly related to the pole 
structure of the S'-matrix. This yields to the expression ( |A.4j ) whose functional form is 
set except for the coefficients a[J' $ appearing in the expansion (|A.8|) of the polynomials 



Qab(@)- The degree fc^fl of these polynomials is fixed by the asymptotic behavior of the 
FF's for large rapidities which depends, of course, on the field $ P[. For the case of 
two-particle FF's of cluster operators, it is easy to see that they are subject to have for 
large 9 at most a constant limit[]. In fact, for two-particle FF's eqs. (|1.6|) read 



Km F%{0) = F? Ff . (2.1) 

->oo 



Hereafter we deal with dimensionless cluster operators which are normalized in such a 
way as to have a vacuum expectation value equal to onef] 

(0|$(0)|0) = F ( f = 1. (2.2) 

In order to fully determine the FF's of the cluster operators we have chosen to focus 
on the set of all one- and two-particle FF's. Listing all the relations among them, one 
obtains a system of equations in the unknown parameters and a^/ $ . Let us see then 
all information we have on the FF's. 

The first equations that one must consider are the dynamical residue equations re- 
sulting from the detailed analysis of the poles they are endowed with. These equations 
relate FF's with different external particles and may have a different origin. In particular, 



4 The limit may vanish in the presence of symmetries. 

5 Since the relevant primary operators will be identified with the cluster ones (except their dimensional 
factors which can be easily restored), in the sequel we will adopt the same normalization also for them. 
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for every simple bound state pole of the amplitude S a b at angle 6 = iu c ah relative to the 
particle A c (see Figure 1), we have 

lim (9 - iu c ab )F*(9) = ir c ab F? , (2.3) 

u— nu c , 

ab 

where the on-mass-shell three-point coupling constant T c ab is given by the residue on the 
pole of the S-matrix 

lim (e- i u c ab )S ab (6)= i (r ab ) 2 . (2.4) 

ab 

Dynamical residue equations are also provided by double order poles and simple order 
poles of type B. Both of them are related to diagrams of the kind shown in Figure 2. For 
each such diagram, one can write the following equation 

lim (9 ab - itfF&e*) = i rSd T%, FfeHl) , (2-5) 

where 7 = 71 — u a cd — u b de . In the case of B poles one can always verify that the amplitude 
S ce (9) has a simple zero at 9 = ij. More complicated residue equations can be in general 
obtained with reference to D poles and higher order ones whose explicit expressions - not 
reported here - can be however easily written, once the corresponding multi-scattering 
diagrams have been identified. 

It must be stressed that the above set of equations just depend on the dynamics of 
the model through its S'-matrix and hold identical for every operator $(a;). Therefore, in 
general, some residual freedom on the parameters is still expected after imposing these 
equations, because they must be satisfied by the FF's of all operators compatible with 
the assumed asymptotic behaviour. 

Adding to this system of linear equations the non-linear cluster equations ( |2.1| ) of the 
two-particle FF's, one obtains in general a redundant set of compatible equations in all 
the unknown parameters of the one- and two-particle FF's. Due to its non-linearity, the 
system allows a multiplicity of solutions which define the so-called cluster operators of 
the theory^. If the number of solutions of the system matches the cardinality of the Kac 
6 In all cases analyzed, the smallest system of equations among different FF's which is sufficient to 
determine their coefficients turns out to involve just a subset of the two-particle FF's. This suggests that 
also in the general case it should be possible to predict the final number of cluster solutions already from 
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table of the model one is led to identify them with the families of FF's of the relevant 
primaries. 

Among the cluster solutions, one can first of all identify the FF's of the deforming 
field <f>i t k- This operator is known to be essentially the trace of the energy-momentum 
tensor O(x) since 

Q{x) = 4vr£ mc lifc , (2.6) 
£ vac being the vacuum energy density which can be easily computed by TBA computa- 



tions 21 



S = (2 7) 

8 EzGPn sm(vrx) 

Here the set P\\ is defined in eq. (|A.l ) and m\ is the lightest particle mass. In view of 
the proportionality fl2.6p , the FF's of <pi^ can be identified among the cluster solutions 
by checking the peculiar equations which characterize the two-particle FF's of G(x) in 
virtue of the conservation of the energy-momentum tensor, namely the normalization of 
the diagonal two-particle FF's 

FZ^) = 2nm 2 a , (2.8) 

and the factorization of the polynomial Q® b for non-diagonal two-particle FF's (a 7^ b) 
into 

Q® fe (coshfl) = (2m a m b cosh 6 + m 2 a + ml) R® b (cosh9) , (2.9) 

where Rf b is a suitable polynomial [I], [J. 

Knowing the FF's of G(x), one is then enabled to make use of the sum-rule ( |1.4|) 
to compute the conformal dimension of the operators defined by the remaining cluster 
solutions in order to identify them with all the relevant primaries of the theory. This 
sum-rule can be evaluated by using the spectral representation of the correlator 

00 „ r\ n f) 

( ©WW ) c = EE / „ 7^^t...,aJ^)^ 1 ,..^(^-0)e-N^ 1 m - osh ^ . 

(2.10) 

In all the models we have studied, the corresponding series for the sum-rule (|1.4j) displays 

a very fast convergence behaviour for any of the cluster operators. The truncated sums 
a "minimal" system, avoiding in this way to deal with systems of equations involving a huge number of 
unknown variables. 
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obtained by including just very few contributions have proved sufficient to attain a good 
approximation of all the values expected by the Kac table of conformal dimensions. In 
this way, the one-to-one correspondence between cluster solutions and primary relevant 
operators can been easily set. 

Finally, having obtained the FF's of all the relevant fields in each integrable defor- 
mation, as a further check of their correct identification, one may employ the formulas 



Ol) relative to the universal ratios of the nearby non-integrable quantum field theories. 



These predictions can then be compared against their numerical estimates obtained from 
Truncated Conformal Space (TCS) approach developed in [I3~|. The agreement between 
numerical estimates and theoretical predictions of the non-integrable effects may provide 
additional confirmation and may remove all possible remaining doubts about the validity 
of the cluster hypothesis for these models. 



3 Integrable Deformations of Ai(2, 7) 

The minimal conformal model -M(2, 7) has, in addition to the identity operator 0^1, only 
two primary operators, 0i 5 2 and 01,3, both of them relevant with conformal weights given 
by —2/7 and —3/7 respectively Q. The perturbations of the conformal action either by 
the "magnetic operator" 12 or by the "thermal operator" 13 are both known to be, 
separately, integrable |Tj|. The S'-matrices and the mass ratios of the two integrable 
models are given in tables Bl and B2. In their massive phase, both perturbations have 
two stable massive particles denoted by A\ and A2, with a mass ratio and a scattering 
matrix which depend on the integrable direction considered. In each case, we expect 
to find two non-trivial independent families of Form Factors solutions to the cluster 
equations ( |1.6| ) (in addition to the family of the null Form Factors relative to the identity 
operator). 

The Form Factors of the primary operators of the model relative to the thermal 
deformation have already been considered in ||. Here, we have performed an ab-initio 
calculation by imposing the cluster equations: our result has been in perfect agreement 
with the FF's of ref. ||, proving in this way that these cluster solutions are also unique. 
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The result of the computation of Form Factors in the two integrable deformations 
[M.{2, 7)] (1,2) and [.M (2, 7)](i i3 ) are summarised in tables 1-2 and 3-4 respectively where 
we list the values of the one-particle FF's and the coefficients a * b \ of the two-particle 
FF's relative to some of the lightest two-particle states. As expected, we find two non- 
trivial solutions of Form Factors families. In each deformation, the FF's of the deforming 
operator suitably rescaled by Q2.6|), can be immediately identified because they satisfy 
the peculiar equations characterizing the trace of the energy-momentum tensor ( p.8[) 
and (|2.9| ). This is further confirmed by employing the spectral representation of the 
correlator (B(x)@(0)) c in the sum-rule ( |1.5| ), which provides in both deformations the 
value of the central charge with a very high precision (the relative error being of order 
10 _4 -10~ 5 ). The identification of both the solutions with the primaries (f>i^ and <pi^ 
is easily established after computing for each solution its UV anomalous dimension by 
means of the sum rule ( |1.4|) . The contributions to this sum rule coming from the dominant 
lightest multiparticle states are given in Tables 5 and 6 for the two deformations (the 
contributions are ordered according to increasing values of the Mandelstam variable s of 
the multi-particle state). The agreement of the truncated sums with the known values 
of the anomalous dimensions is very satisfactory given the fast convergency behaviour 
of the spectral series. In the computation of these sum rules, some three-particle FF 
contributions have been inserted as well, although we do give here their exact expression 
for sake of simplicity (their general parameterization follows the one adopted for instance, 
in [Q]). It should be noticed that the oscillating behaviour of these sums is typical of non- 
unitary theories where one expects, in general, both positive and negative terms. 



3.1 Non— Integrable Deformations of ,M(2,7) 

For each possible integrable deformation of the model, the addition of a further orthog- 
onal deformation breaks its integrability leading, among other things, to corrections of 
the mass spectrum and of the vacuum energy. Both corrections can be independently 
computed by performing a numerical diagonalization of the off-critical Hamiltonian by 



means of the so-called Truncation Method ||13|| . We have carried out this analysis com- 
paring these non-integrable data with the theoretical predictions by eqs. (|1.3|). Let us 
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briefly describe the output of these studies. 
The double non-integrable deformation 

[M(2,7)] (lj3 ) + £0i )2 , 

for small values of emf 1,2 2 has already been studied in f9j, where a good agreement 
between numerical and theoretical values has been found. Having obtained the FF's for 
the 0i 5 2 deformation, we are now able to complete the analysis by testing the opposite 
deformation 

[^(2,7)1(1,2) +e0i, 3 • 

The numerical determination of the two universal ratios of eq. (|1.3|) (for small values 
of emf 1,3 2 ) gives = 0.675 and = — 0.244 mf with a precision estimated to 
be up to a few percents. This values fully agree with the computed theoretical values 
= 0.68404 and = -0.24365 mf' (see, for instance Figure 4 and 5 where the data 
relative to the ratios I 221 and +^ respectively are reported for different values of e). 

om2 omi / 

4 Integrable Deformations of A4(2,9) 

In this section, we turn our attention to the Ai(2,9) minimal model which displays a 
richer structure in the RG space of relevant couplings. This model has in fact, besides the 
identity, three primary operators (fti^, <j>\ z an d </>i,4 which are all relevant with conformal 
dimensions —1/3, —5/9 and —2/3 respectively. These fields taken separately give rise to 
different integrable deformations of the conformal model, each of them characterized by 
a different mass spectrum and ^-matrix (see tables B3, B4 and B5 in Appendix B). In 
particular, the first two deformations produce three-particle mass spectra (with different 
mass ratios) while the last one gives a four-particle spectrum. 

The FF's of the primary operators in the ^^-deformation had already been obtained 
in ref. || and were known to satisfy the cluster property. Again, our derivation of these 
FF's as solutions of the cluster equations proves that the FF's found in || are the only 
possible cluster solutions. 

The Form Factors of the cluster solutions for each of the three above mentioned 
deformations have been computed according to the strategy explained in Section 2. The 
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resulting one-particle FF's and two-particle FF's coefficients are given in tables 7-8, 9-10 
and 11-12 respectively. The important result is that in each integrable deformation of 
this model, three families of non-trivial solutions have been found. Among the solutions, 
we have firstly identified the FF's of the deforming field by checking the exact fulfillment 
of eqs. ( |2.8|) and ( |2.9|) , after the appropriate rescaling (|2.6|) . Moreover, the c-sum-rule 



OP can be easily shown to give very precise approximations of the central charge in 
each of the three separate deformations. 

As for the other solutions, they have been successfully identified with the FF's of 
the primary operators by computing their anomalous dimension by means of eq. ( |l-4j ). 
The first contributions to these sums are given in tables 13, 14 and 15. In all cases the 
agreement with the expected anomalous dimensions of the primaries is established, even 
though the convergence of the series is observed to be noticeably faster for lower absolute 
values of the anomalous dimension of the deforming field. This observed trend is indeed 
expected from the short-distance behavior of the correlator (|2.10|) , as predicted by the 
Operator Product Expansion of the fields. In fact, in the models A / l(2, 2 n + 1) where the 
fields have negative anomalous dimensions, this correlator displays a zero at the origin 
whose power law exponent is larger for lower absolute values of the anomalous dimension 
of Q(x); correspondingly, the small x region of integration in ( |1.4| ) is less relevant making 
the lightest multiparticle states more dominant in the series. 

4.1 Non— Integrable Deformations of A4(2,9) 

The availability of the FF's of all the primary fields of the model has allowed us, in 
each of the three separate integrable deformations, to consider two different orthogonal 
non-integrable deformations. We have had then the possibility of testing the theoretical 



values obtained for the universal quantities ( |1.3|) versus their numerical TCS estimates 
in six different multiple deformations, exploring in this way the non-integrable region 
around the conformal point of the model. The outcome of the analysis in all the defor- 
mations is summarized in Table 16. Since the precision of TCS data is expected to be of 
approximately a few percents, the comparison with the computed theoretical values is in 
all cases quite satisfactory. 
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5 Conclusions 



The main purpose of this work has been to substantiate by means of concrete ab-initio 
calculations the cluster hypothesis for the Form Factors of the relevant operators in in- 
tegrable quantum field theories obtained as deformation of a conformal action. We have 
studied, in particular, the matrix elements of the primary operators in the integrable 
deformations of the first models of the non-unitary series Ai(2,2n +1). In all cases 
analysed, we have confirmed the cluster hypothesis since we have found a one-to-one 
correspondence between the independent solutions of the cluster equations and the rele- 
vant fields. 

It should be said that the absence of internal symmetries of the above models has 
played an important role in carrying out our computations. In fact, in this situation one 
can exploit the cluster equations Ql.Ej ) in their full generality. It would be interesting to see 



how the results of ref. |T0| generalize to the case of quantum field theories with internal 
symmetries which induce selection rules on the matrix elements. Another important 
open problem is also to understand the meaning of the cluster properties in quantum 
field theories which cannot be regarded as deformation of conformal models. A complete 
understanding of all these aspects of the Form Factors would allow us to better understand 
the asymptotic high-energy regime of quantum theories and their operator content. 

Acknowledgements. We are grateful to G. Delfino and P. Simonetti for useful discus- 
sions. 
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Appendix A 



In this appendix we give the general parameterization adopted throughout the paper for 
two-particle Form Factors. The ^-matrices of the specific models analysed in this paper 
are given in Appendix B where the generic amplitude 

s*(e)= n (*) Px II (~y) qy (a- 1 ) 

xeP ab yez ab 

is written adopting the notation 

. , tanh \ (9 + ma) . . . 

(a) = -f-^ : J r A.2 

v ; tanh \ (9 - tira) 

and the positive rational indices x and y label the poles and the zeros displayed by the 
amplitude in the physical strip Im^ G [0, 7r]. The bound state simple poles of the S- 
matrices are identified by superscripts which denote the particles produced. There are 
however also simple poles of a different nature which have been labelled with suffices B 
and T> that are not related to any bound state and are due to multiparticle scattering 
processes of the kind shown in Figure 2 and 3 respectively. The fact that these diagrams 
(which usually produce second and third order poles) are here responsible for simple 
poles is due to the occurrence of zeros in the S"-matrix factors carried by the internal 
crossing lines [ 20| . Higher order poles are present as well and, among these, in the model 



[«M (2/9)] (1,4), some triple poles labelled by an asterisk which also have a non-standard 
diagrammatic interpretation. The understanding of the nature of all the poles is necessary 
in order to assign to the FF's the correct pole structure. 

The general two-particle Form Factor of a scalar operator $(x) 

*2 fa " *2) = ( | $(0) | AM A b (9 2 ) ) , (A.3) 
will be parameterized by 



L>ab (?) 



where the "minimal" FF 



n«r«o 

F„f (0) = (-isinh W 2)) » |f — , (A.5) 



y£Z ab 
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which has neither zeros nor poles in the physical strip, is written in terms of the function 

fc+i 
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(A.6) 



This function is normalized by g x (in) = 1 and behaves asymptotically as 

g x {6) ~ e |e|/2 for -> oo . 



(A.7) 



The factor Q® b (6) in ([A. 4 ) is a polynomial in cosh carrying the dependence on the 
specific operator 

Qa,= £a«,COSh*(0). ( A - 8 ) 
k=0 

The most subtle element in the parameterization of the FF's is represented by the struc- 
ture of the poles which, in eq. ( |A.4| ) are introduced by the factor D ab {6). In order to 
establish which poles are to be found in a FF one must in general have a complete un- 
derstanding of the nature of the poles in the corresponding S-matrix element in terms 
of microscopical processes. We will write in general, 



D ab (9) = n K x vt, 



(A.9) 



where the set of indices is defined in ( |A.1|) and 

cos(7nr) — cosh(6 l ) 



V x (9) 



2cos 2 (f ) 



(A.10) 



For bound-state simple poles and ordinary higher order poles of the S— matrix, the correct 
rule for determining the indices i x and j x is given by H 



x Tl , J x 
X 1>> ) Jx 



(All) 



n — 1 , if p x = 2n — 1 ; 
n , if p x = 2n , 

For simple poles of type (x)g and (x)u one can show that the correct indices are still 
i x = 1 and j x = 0, as for a bound state simple pole. 

Notice however that the poles of the FF's induced by the triple poles labelled with 
* in [M. (2/9)] (1,4) do not fall within the above analysis. Their general expressions is not 
further investigated here since these FF's were not needed in the present work. 
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As a final remark, notice that every function (a) could be equivalently written as 
(1 — a) without changing the S'-matrices. However, the pole prescription given above 
for the FF's is sensitive to this change in the case of odd order poles. Therefore, all the 
labels a in the ^-matrices reported here have been chosen to give (in units of in) the 
value of the direct s-channel resonant angles in the case of bound state odd poles and 
also in the case of poles of type B and D R. Only with this choice, the above prescription 
gives the correct poles of the FF's. 

With the parameterization ( |A.4| ), the two-particle Form Factor of a general operator 



$ is therefore completely determined after fixing the coefficients \ in the expansion 



7 For s-channel in these cases we mean the one defined by Figures 2 and 3 with particles flowing 
upward 
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Appendix B 



In this Appendix we give the S'-matrices of the integrable models analyzed in this work. 
The function (a) used in the tables is given in eq. ( |A.2[ ). Anomalous simple poles have 
been labelled with B and V, while the anomalous triple poles of the model [A / i(2/9)](i j4 ) 
are identified with *. 



Sn(0) = (1) (|) (-1) 
Sn(f» = (§) (S), 

s-W = (I) (I) B (I) 



m 2 = 2 cos — mi = 1.9696... mi 



Table Bl: S-Matrix and mass ratios of the [M(2/7)]n ) 2) model. 



Sn(9) = (I) 

= (!) (i: 


to 


71 

rri2 = 2 cos — mi = 1.61803 . . . mi 

5 



Table B2: S-Matrix and mass ratios of the [M(2/7)]n i 3) model. 
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Q I fl\ — 

&ii\y) — 


Am / i\ 


£>12\V) — 


1 3 

f 23\ f l\ f 5\ ( 5 \ 
V 24 ^ V 8 / WgV 24 y 


>J13\V) — 


\12) \12) B 


C* { Q\ 
022{t>) = 


2 „ 
/2\ / 11 \ 2 ( 7\ ( l\ 
\3j \12J \l2j v V 4j 


S23(@) = 


(i) (1) (1) (i) 

\24/ \8/b Vo/g \24Jx> 


5*33 (#) = 


(§) (12)0 (l) B (n)x, (2)2? 




sin^ 

Hfrni a = l, 2, 3 
smfj 


m a = 



Table B3: S-Matrix and mass ratios of the [M(2/9)](i j2 ) model. 



Su(6) = 


2 

(?) 




Sn(0) = 


(?) 


3 

(f) 


Sis(9) = 


2 

W 


3 

f 4\ 


S22{9) = 


3 

(?) 


(I) 2 


S23(9) = 


(?) 


(I) (f) 2 


SW^) = 


(?) 


(f) 2 (f) 2 


sin^ 

m a = . ; 

sm y 


mi 


a = 1, 2, 3 



Table B4: S-Matrix and mass ratios of the [M(2/9)](i j3 ) model. 
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On 


\ ) — 


(l\(7\(l W 1 W 2\ 
UJ UJ UJ I 15 J I 5 J 


5*12 


[9) = 


1 3 

( 23 \ / 13 \ 
\30j \3o) 


5*13 


'Q\ — 


(^) (?A (I) f3\ 2 / 2 W l\ 






fi^ m (i\ 2 

\15j B \3j 


5*22 


.17 J — 


(A (A ( 8 ) 
\3j \5J \l5j B 


Q 

023 


— 


{6) \2) B \\0) B \30j B 


024 


\ ) — 


(9\ (23\ (3\ (19\ (17\ 2 
\10) \30j B \10) B \30j B \30j 


033 


v^) — 


3 

/2\ 3 ( 2 \ 2 ( 7 \ 2 ( 1 "\ / 2\ 
(3) (is) (ij (-15) (-5) 


5*34 


[9) = 


(A) ft) (H)* (if 


S44 


[0) = 


4 

(3)* (is)* (s)x> (15)0 (is)b (s) 




m 2 


7tt 

= 2 cos — mi = 1.48629... mi 
30 




m 3 


7T 

= 2 cos — m x = 1.95630... mi 
15 




m 4 


= 2 cos ^ m 2 = 2.82709... mi 



Table B5: S-Matrix and mass ratios of the [M(2/9)](i j4 ) model. 
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Table Captions 



Table 1 One-particle form factors of cluster solutions in [_M(2/7)](i )2 ). 

Table 2 Two-particle form factors coefficients of cluster solutions in [A / f(2/7)](i j2 ). 

Table 3 One-particle form factors of cluster solutions in [M (2/7)] (i j3 ). 

Table 4 Two-particle form factors coefficients of cluster solutions in [A / f(2/7)](i i3 ). 

Table 5 Sum rules of the conformal dimensions of primary operators in [A / f(2/7)](i )2 ). 

Table 6 Sum rules of the conformal dimensions of primary operators in [A^(2/7)](i j2 ). 

Table 7 One-particle form factors of cluster solutions in [_M(2/9)](i j2 ). 

Table 8 Two-particle form factors coefficients of cluster solutions in [_M(2/9)](i j2 ). 

Table 9 One-particle form factors of cluster solutions in [M.(2/9)]^ lj3 y 

Table 10 Two-particle form factors coefficients of cluster solutions in [M(2/9)]( lj3 y 

Table 11 One-particle form factors of cluster solutions in [M. (2/9)] (1,4). 

Table 12 Two-particle form factors coefficients of cluster solutions in [_M(2/9)](i j4 ). 

Table 13 Sum rules of the conformal dimensions of primary operators in [A / f(2/9)](i ;2 ). 

Table 14 Sum rules of the conformal dimensions of primary operators in [A / f(2/9)](i ;3 ). 

Table 15 Sum rules of the conformal dimensions of primary operators in [_M(2/9)](i j4 ). 

Table 16 Comparison between numerical and theoretical estimates of data obtained in 
different non-integrable deformations of -M( 2) g). . 
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o 


01,2 


01,3 


F° 


0.8129447456 i 


1.245503611 i 


F° 


-0.1200387686 


-0.4656766285 


Table 1 


O 


01,2 


01,3 


a u,o 
a u,o 


-0.6905355776 
1.570496171 


-0.4178217785 
3.686419944 


a 12,0 

a {1) 

u 12,0 


31.91217166^ 
25.76337182 i 


160.8200658 i 
153. 1262244 i 


a (0) 

u 22,0 
"22,0 


-12.74804909 
-3.97663589 


-71.66459155 
-59.84689851 


Table 2 


O 


01,2 


01,3 


F° 


0.8703387193 i 


1.408237641 i 


F 2 ° 


-0.3322661173 


-0.8698840033 


Table 3 


O 


01,2 


01,3 


r7 (0) 
"11,0 


-1.453085043 


-3.804226098 


a (0) 
"12,0 

a (1) 


10.38924846 i 
6.420908640 i 


30.40986050 i 
27.19940617 i 


a (0) 

a 22,C 
(1) 

a 2 2,o 


-13.76381909 
-4.702281947 


-42.18951412 
-32.22992104 



Table 4 
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states 


s 


A 12 -terms 


A 13 -terms 


A 1 


1.000 m! 


-0.2922910 


-0.4478157 


A 2 


1.969 mi 


0.0016428 


0.0063729 


A! A, 


> 2.000 mi 


0.0051123 


0.0137590 


A 1 A 2 


> 2.969 mi 


-0.0000763 


-0.0004400 


At A x Ai 


> 3.000 mi 


-0.0001040 


-0.0004777 


A 2 A 2 


> 3.939 mi 


0.0000003 


0.0000040 


sum 


-0.2857159 


-0.4285976 


value 


expected 


-0.2857143 


-0.4285714 



Table 5 



states 


s 


A12 -terms 


A 13 -terms 


A 1 


1.000 mi 


-0.3221795 


-0.5212974 


A 2 


1.618mi 


0.0290206 


0.0759768 


A l A 1 


2.000 mi 


0.0098699 


0.0258398 


A l A 2 


2.618 mi 


-0.0023149 


-0.0089996 


A 1 A 1 A 1 


3.000 mi 


-0.0003334 


-0.0013803 


A 2 A 2 


3.236 mi 


0.0001155 


0.0006612 


sum 


-0.2858218 


-0.4291998 


value expected 


-0.2857143 


-0.4285714 



Table 6 



o 


01,2 


01,3 


01,4 


F° 


0.7548301717 i 


1.288575652 i 


1.564862744 i 


if 


-0.1056909725 


-0.4593398099 


-0.7331609072 


if 


-0.01375684037 i 


-0.1 175389994 i 


-0.2854817817 i 



Table 7 
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o 


01,2 


01,3 


01,4 


a {0) 
a n,o 

a u,o 


-0.3810248990 
1.289925788 


0.1280888115 
3.759118917 


0.6449629545 
5.543942595 


a (0) 

a 12,0 

a {2) 
a i2,e> 


14.10905183 i 
-12.74323779 i 
-19.36998044 i 


75.18632019 i 
-79.90895489 i 
-143.7096872 i 


110.3472056 i 
-180.9845092 i 
-278.5592522 i 


a 13,C> 
a 13,C 


-1.826322080 
-1.116015559 


-18.97540047 
-16.27774386 


-51.56786333 
-48.01279071 


a (0) 

"22,0 
a 22,C 


-1.466545085 
7.821352950 
2.717967823 


-3.003367424 
60.49540624 
51.33773403 


14.9160654 
160.4007705 
130.7877664 


a 23,C 

a 23,C 
(2) 

a 23,e> 


153.8279467 i 
175.5584268 i 
30.43124786 i 


1842.946063 i 
2962.508857 i 
1130.002086 i 


5426.66381 i 
9796.436391 i 
4380.673323 i 


a 33,0 

a (1) 
tt 33,0 


-32.42110324 
-20.23293766 
-2.174915595 


-450.0936155 
-589.1376530 
-158.7701993 


-1394.808207 
-2309.626757 
-936.6165096 



Table 8 






01,2 


01,3 


01,4 


F° 


0.8020765716 i 


1.445292066 i 


1.802249672 i 


Ff 


-0.3139111339 


-1.019263084 


-1.584911324 


F 3 ° 


-0. 1373692453 % 


-0.5561967434 i 


-1.002231818 i 



Table 9 
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o 


01,2 


01,3 


01,4 


a (0) 
a n,e> 


-0.9631492344 


-3.127326026 


-4.862860736 


a 12,0 


10.64696613 i 


40.73951464 i 


72.35568181 i 


a (1) 

u 12,0 


5.908620424 i 


34.57048356 i 


67.03219861 i 


a 13,C 


-2.592348236 


-11.32977918 


-21.24912975 


a (1) 
u i3,e> 


-1.153703500 


-8.417302355 


-18.91350458 


a (0) 


-5.978990567 


-26.44069921 


-49.87674876 


a 22,0 


-1.544771430 


-16.28633559 


-39.37864116 



Table 10 



o 


01,2 


01,3 


01,4 


Ff 


-0.9043544898 i 


-1.72785339 % 


-2.211259663 i 


Ff 


-0.5483648961 


-1.476188315 


-2.169493373 


Ff 


0.2673316508 i 


0.8709319528 i 


1.45902371 i 


F? 


-0.08488118964 


-0.3489749771 


-0.6451795597 



Table 11 
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o 


01,2 


01,3 


01,4 


a {0) 
a n,o 

a n,o 

(2) 

a n,o 


1.623982681 
-0.9778411563 
-2.029027259 


4.256426530 
-1.325966569 
-7.406691607 


6.219867507 
-1.477684504 
-12.13081476 


a 12,0 
u 12,0 


-9.935037127 i 
-4.790105254 i 


-30.48935000 % 
-24.63686052i 


-50.02403946 i 
-46.33773309 i 


a 13,C> 

a 13,0 

(2) 
a 13,C 

(3) 
a 13,0 


-45.21074197 
-441.3756086 
-533.3237140 
-139.4173540 


-145.9600730 
-1583.189947 
-2301.408527 
-867.7984505 


-220.9609710 
-2731.357697 
-4364.089257 
-1860.500753 


a 14,0 

a 14,0 

(2) 
a 14,0 


44.34032961 % 
54.79194008 i 
11.43395882 i 


189.6077348 i 
275.9507675 i 
89.81474554 i 


357.8357762 i 
562.9638801 i 
212.5038539 i 


(0) 

a 22,e> 

"22,0 


-9.190266093 
-2.709639668 


-30.91094092 
-19.63612453 


-52.20534546 
-42.41201502 


(0) 

a 23,e> 

a 23,0 

fl (2) 
"23,0 


-81.75802420 i 
-92.61128143 i 
-16.66533965 i 


-304.2244838 i 
-446.8601169 i 
-146.1571720 i 


-530.8872409 i 
-884.9723034 i 
-359.8444599 i 



Table 12 
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states 


s 


A12 -terms 


A13 -terms 


A 14 -terms 


A, 


1.000 mi 


-0.3409847 


-0.5820972 


-0.7069063 


A 2 


1.982 mi 


0.0017003 


0.0073894 


0.0117945 


A 1 A ± 


> 2.000 mi 


0.0061957 


0.0207909 


0.0316698 


A 3 


2.931 mi 


-0.0000132 


-0.0001126 


-0.0002734 


Ai A 2 


> 2.982 mi 


-0.0000951 


-0.0007084 


-0.0014392 


A x A x A 1 


> 3.000 mi 


-0.0001421 


-0.0009038 


-0.0017386 


A 1 A 3 


> 3.931 mi 


0.0000009 


0.0000117 


0.0000339 


A 2 A 2 


> 3.965 mi 


0.0000004 


0.0000061 


0.0000157 


A 2 A 3 


> 4.914 mi 


-0.0000000 


-0.0000002 


-0.0000008 


sum 


-0.3333379 


-0.5556241 


-0.6668445 


value expected 


-0.3333333 


-0.5555556 


-0.6666667 



Table 13 



states 


s 


A 12 -terms 


Ais-terms 


A i4 -terms 


A x 


1.000 mi 


-0.370679 


-0.667941 


-0.832909 


A 2 


1.802 mi 


0.031509 


0.102310 


0.159088 


A, A, 


> 2.000 mi 


0.013898 


0.045127 


0.070170 


A 3 


2.247 mi 


-0.004839 


-0.019592 


-0.035304 


At A 2 


> 2.802 mi 


-0.003604 


-0.018722 


-0.035573 


A x A 1 A x 


> 3.000 mi 


-0.000628 


-0.003514 


-0.006763 


A x As 


> 3.247 mi 


0.000663 


0.004114 


0.008844 


A 2 A 2 


> 3.604 mi 


0.000211 


0.001684 


0.003864 


sum 


-0.333469 


-0.556534 


-0.668583 


value expected 


-0.333333 


-0.555556 


-0.666667 



Table 14 
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states 


s 


A 12 -terms 


A 13 -terms 


A14 -terms 


A 1 


1.000 mi 


-0.451081 


-0.861833 


-1.102950 


A 2 


1.486 mi 


0.121478 


0.327017 


0.480603 


A 3 


1.956 mi 


-0.022989 


-0.074895 


-0.125468 


A l A 1 


> 2.000 mi 


0.035896 


0.121577 


0.197637 


Ai A 2 


> 2.486 mi 


-0.023279 


-0.101138 


-0.183618 


A, 


2.827 mi 


0.001546 


0.006354 


0.011748 


A 1 A 3 


> 2.956 mi 


0.004304 


0.022374 


0.045474 


A 2 A 2 


> 2.973 mi 


-0.001535 


-0.009929 


-0.022429 


A 2 A 3 


> 3.443 mi 


-0.000330 


-0.002101 


-0.004686 


A 1 A 4 


> 3.827 mi 


0.003595 


0.020054 


0.040870 


sum 


-0.332396 


-0.552519 


-0.662819 


value expected 


-0.333333 


-0.555556 


-0.666667 



Table 15 





Srrii 
5m 2 




deformation 


numerical (±3%) 


theoretical 


numerical (±3%) 


theoretical 


[7W(2,9)] ( i, 2 )+£0i,3 


0.590 


0.592049 


-0.275 


-0.275404 


[M(2, 9)1(1,2) +£0 M 


0.661 


0.660963 


-0.204 


-0.204124 


[A<(2, 9)1(1,3) +e&,2 


0.390 


0.391396 


-1.04 


-1.03826 


[M(2, 9)](i,3) +£0 M 


0.811 


0.83681 


-0.205 


-0.205640 


[M(2,9)] (1 , 4 )+£0l,2 


-0.133 


-0.131367 


1.73 


1.74582 


[A^(2,9)] {1 ,4)+£01,3 


0.238 


0.240486 


-0.550 


-0.548156 



Table 16 
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Figure Captions 



Figure 1 Bound-state simple pole diagram. 
Figure 2 "Butterfly" diagram. 
Figure 3 "Dragonfly" diagram. 

Figure 4 Numerical TCS estimates of 8mi versus Sm2 for different values of the "non- 
integrable" coupling in the model [_M(2/7)](i )2 ) +£0i,3- The dashed line represents 
the theoretical prediction. 

Figure 5 Numerical TCS estimates of SS vac versus mf^brrii for different values of the 
"non-integrable" coupling in the model [A4(2/7)](i ;2 ) + £0i,3- The dashed line 
represents the theoretical prediction. 
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Figure 4 
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